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Abstract 

We consider a multi-species generalization of the Asymmetric Simple Exclusion Process on an open chain, 
in which particles hop with their characteristic hopping rates and fast particles can overtake slow ones. 
The number of species is arbitrary and the hopping rates can be selected from a discrete or continuous 
distribution. We determine exactly the phase structure of this model and show how the phase diagram 
of the 1-species ASEP is modified. Depending on the distribution of hopping rates, the system can exist 
in a three-phase regime or a two-phase regime. In the three- phase regime the phase structure is almost 
the same as in the one species case, that is, there are the low density, the high density and the maximal 
current phases, while in the two-phase regime there is no high-density phase. 
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1 Introduction 

The asymmetric simple exclusion process ^, || ^ ||, ^ refers to a collection of Brownian particles which 
under the influence of a driving force, do biased random hopping on a one dimensional lattice and interact 
via hard core repulsion with each other. In the totally asymmetric case each particle is injected to the system 
from the left with rate a and hops only to the right neighboring site with a rate normalized to unity and 
finally is extracted at the right end with rate /3. 

This is a model far from equilibrium with many re-interpretations which makes it a suitable model for 
studying such phenomena as diverse as sur face growth 0, and traffic flow g |], (see [|o|,|l| |T| and 
references therein.) 

One of the most interesting aspects of this process is the possibility of boundary induced phase transitions. 
It has been observed through various types of solutions |l^, |l^, |l^ that by changing the rate of 

injection and extraction of particles, different phases will develop in the system. The phase diagram of the 
model representing the macroscopic current of the particles in various domains in the a — (3 plane is depicted 
in fig.(I). 

Although some of the characteristics of the process have been obtained by other analytical methods,]^ 
[0,||l8|, the solution by Matrix Product Ansatz (MPA) has proved much useful for obtaining among 
other things, the same results in much simpler ways. The general formulation of MPA [|T9| has also been 
shown to be amenable to further generalizations ||2^, U, |2^, |2^, |2^ . 

In the MPA the probabilities P(ti • ■ -tl) , where is the random variable associated with site i (being 
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Figure 1: Phase diagram for the single species ASEP 



for empty site and 1 for an occupied site), is written as: 



1 ^ 

Pin, ■■■tl) = — <W\ WinD + (1 - n)E\V > 



where the operators D and and the vectors < W\ and \y > satisfy the following relations: 



(1) 



DE=D+E 



(2) 



D\V>^-\V> 



< W\E =- <W\ 
a 



(3) 
(4) 



In (|^), Zl is a normalization constant and is suitably called the partition function. Its value is given by 



Zl =< W\C^\V > 



(5) 



where C := D + E. 



The question of a natural species generalization of the ASEP, so that in the special casep = 1, one obtains 
the results of the one species ASEP has been answered in the affirmative in |2^, by postulating 
a generalization of the algebra (Q-^), which we call the p— ASEP algebra. Since the number of species is 
quite arbitrary and it can even be infinite, in which case the hopping rates are taken from a continuous 
distribution, we rewrite the multi-species algebra in a more general form than that of [ p2[ . This new algebra 
is generated by a discrete generator E and a one parameter family of generators D{v) where v £ These 
generators are subject to the following relations: 

D(v)E = -D(v)+E (6) 

V 

n( '\nf ^ vD{v')-v'D{v) , 

Div )0(v) — ; V > V (7) 

V — v' 

D{v)\V> = J \V> (8) 
w + p - 1 

< W\E = - <W\ (9) 
a 

We also need to show how this continuously parameterized algebra is derived in the MPA formalism. This is 
done in the appendix. The hopping rates v are taken from a general distribution a{v), with support [ui, oo). 
That is, Vi is the smallest hopping rate in the ensemble. Note that this ensemble refers to the particles 
waiting to enter the system, or the ensemble of particles moving in the system, if there were no interactions 
between the particles. For this reason we call v the intrinsic hopping rate or average velocity of a particle. 
The process described by this algebra is one in which each particle of velocity v arrives at the left end with 
rate a(v) := av, (i.e: the input current is (T{v)av), hops to its right neighboring empty site with rate v and 
leaves the system at the right end with rate l3{v) — P + v — 1. If this particle encounters on its way a site 
occupied by a particle of intrinsic velocity v' with v' < v, it will overtake it with rate v — v' , otherwise it 
stops. For all the extraction rates f]{v) to be positive we also require that (3 > 1 — vi. The unit of time is set 



so that the average hopping rate is unity, i.e: J va{v)dv = 1. Thus the parameters a and /3 are respectively 
the total injection rate and the average extraction rate of the particles respectively. Note that although 
the multiplicative dependence of the injection rate a{v) on v is rather natural, this is not the case for the 
extraction rate f3{v) the form of which is dictated only by our demand to solve the system through the MPA 
relations (|[^)- Note also that all the elementary processes are stochastic, i.e: in a time interval dt, a particle 
of velocity v present in a given site, hops to the right empty site with probability vdt, and does not move 
with probability 1 — vdt. 

The model we consider depends on two boundary parameters a and (3 and on the distribution function a{v). 
All our arguments below are also valid for a discrete distribution, for which a{v) i J^i ^("^ ~ '^O- 

The main motivation for pursuing this problem is to see how the phase structure of the one-species ASEP 
(hereafter denoted by 1-ASEP) is modified, when we have particles with a variety of hopping rates and es- 
pecially when particles can overtake each other. Do we still have the phases of low-density, high-density and 
maximum Currents, present in the 1-ASEP, or is it changed in an essential way?. How does the variety of 
hopping rates in the bulk or their probability distribution enter the picture and what role does this distri- 
bution play in the phase structure of the system ? How the absence of particle-hole symmetry in this model 
is reflected in the phase diagram? 

We will go through these questions by providing an exact solution of this problem, and will obtain a gener- 
alization of the phase diagram of the 1-ASEP. 

As far as we consider only the mean field line a + /3 — 1 , one dimensional representations of the algebra 
(^[^ give an exact solution. However to uncover the important role of fluctuations, we should explore the 
full a — f3 plane and for this we should use the infinite dimensional representation. What we will do is 
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to calculate exactly the generating function for partition functions of systems with different sizes and by 
carrying out an analysis of its singularities, determine the currents and the different phases of the system. 
The phase structure depends on the values of a, /3, and on the characteristics of the distribution function. 

Main Results: 

• To every distribution function a{v) of hopping rates, we can assign a real number l[a], defined in (^7|), 
which essentially depends on the behavior of a{v) for small hopping rates (i.e: if cr(wi) = or not and if yes 
how slowly it approaches this value). The parameter Ic — is special in the sense that for all distribution 
functions with l[a] < 0, the phase diagram of the multi-species ASEP is almost the same as the phase dia- 
gram of 1-ASEP, that is, in the a — f3 plane we have three regions of low density, high density and maximum 
current phases. The value of the maximum current and the shape of the coexistence curves between different 
phases depend on the distribution function (see figs.(2)and (3)). We also obtain the average density of all 
types of particles in all three phases. 

• If on the other hand l[(j] > 0, then the phase diagram consists of only two phases, namely the low 
density and the maximum current phase. The extraction rate f3 does not have any effect on the system 
and only the injection rate a determines which phase will develop in the system (fig. 4). We also obtain the 
average density of all types of particles in both phases. 

Thus the general shape of the phase diagram is controlled by three parameters. The two control pa- 
rameters a and (3 represent the effect of boundaries and the third parameter Z[cr], takes into account the 
distribution function of hopping rates. 
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Figure 2: Phase diagram for multi-species ASEP when l[a] < 



This paper is organized as follows: In section (2), wc introduce the necessary preliminary material from p2[ 
and set our notations and conventions. In section (3) the generating functions for the currents and average 
density of particles of each species are introduced and the former one is calculated exactly. In section (4) the 
total current is calculated in terms of which the different phases of the system are determined. Section (5) 
is devoted to the calculation of the generating function for the average densities and the calculation of latter 
quantities. In section (6) we discuss two special cases, namely the 1-ASEP where we reproduce the already 
known results, and the case when one of the hopping rates is much smaller than the others. We conclude 
the paper with a discussion in section (7) . 
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Figure 4: Phase diagram of the multi-species ASEP when /[cr] > 



2 Some algebraic preliminaries 



To make the present paper self-contained, we quote the basic definitions and theorems from [ p2[ , to which 
the reader can refer to, for further details and proofs. 



The algebra (jg-^ has only one or infinite dimensional representations. We write the infinite dimensional 
representation in a convenient basis consisting of vectors |0 >, |1 >, |2 >, ■ • ■, with the following actions: 



E\n > = In + 1 > 



D{v)\n> = v-''\0>+v-''+'\l>+v-"+^\2>+...v-'\n-l>+\n> 



(10) 
(11) 



We also have: 



|F>=£(f-/3nn> 



(12) 



rs=0 



<W\ = Y^ia)'^ < n\ 

n=0 

The basis{< n\} is the dual of the basis {\n >}, i.e: < 7i|?7i >= Sn,: 



(13) 



The operator C is defined as 



C ■.^E+ / aiv)D{v)dv 



(14) 



whose action on the basis vectors is calculated to be 



n+l 



C|n ^ < -1^ > Ifc > 



(15) 



k=0 



where < > is the average of the inverse fc— th power of the hopping rates. (< >= ^ Ej=i (v])'' ^'^^ 
discrete distributions and < >= / dva{v)-i^ for continuous distributions). The explicit matrix form of C 



:s : 



C 



1 < r: > < TT? > < 73 > <-ki> 



1 1 <^> <^> <^> 



\ 



1 











1 < i > 



(16) 



Clearly direct evaluation of the N-the power of this matrix is a formidable task. In the one species case 
where < >— 1, Vfc, one can go to a new basis {\n >'.— {E — 1)"|0 >} in which the matrix C becomes 
tri-diagonal with the simple form Cn,m — '2Sn,m + Sn,m+i + Sn,m-i- This Hermitian matrix can then be 
either easily diagonalized |0| or else, yields simple recursion relations which can be solved by an analogy 



with the master equation of a random walk in the presence of an absorbing wall [l^ . Due to the complicated 



form of the matrix C, none of the above strategies work in the present case. The same is true with the 
method of repeated appUcation of the algebraic relations (||-^ and calculating directly the matrix element 

< w\c'^\v > 0. 

There is however one basis in which a manageable recursion relation can be found, namely the coherent basis 
defined as follows: 

oo 

|U>=^U"|71> (17) 
n=0 
oo 

<u|==^w"<n|. (18) 

n=0 

Note a slight difference of our notation with that of |^ in the symbol for < u|, where this state would have 
been denoted by < u~^\. These states have the following properties: 

<u\E^u<u\ (19) 

D(v)\u >^ —^\u > (20) 

V ~ u 

< uj\u >= — - — for \uuj\ < 1 (21) 

1 — ULO 

\u><u-'\^I (22) 



27rm 

where / is the identity operator and c is any contour encircling the origin. 

Remark: In calculation of matrix elements of operators between two states < a\ and \b >, one can insert 
any numbers of unit operators in the form of ( |2^ ) with integration variables ui,U2, ■ ■ ■ from left to right, 
provided that |i| > > \u2\ > ■ ■ ■ \b\. This is due to the restriction (|2l|). The results of such calculations 
are then valid only for \ab\ < 1 and must be analytically continued to larger domains. From the definition 
of C and (|9l)-(|l|) one obtains 



< u\C\w >=(u+ [ dvaiv) ^ ) (-^) =: U + 9{w)) {-^—) 



(23) 
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where the second equahty defines the function g{'w), a shorthand and useful notation for which is 



Here the average is taken with respect to the probabihty distribution of hopping rates, (j{v). 



3 The generating functions for current and average densities 

The total current for a system consisting of N sites has been found to be ||2^ 

<a\C^~^\b> 

~ <a\C^\h> ' ^^^^ 
where < a\ and \h > are coherent states and for convenience, we have denoted 1 — /3 by 6, and by a. 
The current density of particles of velocity v is given by 

J{v) = va{v)J (26) 



In the thermodynamic limit — > oo, there is a simple way to evaluate (25). Define a generating function 



/(s;a,6) ^ < a|C^|5> (27) 



The convergence radius of this formal series, R, is precisely what we need. In fact 

w^oo <a|C^|&> 

Remark: The function /(s; a, b) has a Taylor-series expansion in terms of (non-negative) powers of its three 
arguments, which means that there is a region containing the origin of the space C"^, where / is analytic. 
This is easily seen by noting that can in principle be expanded in terms of integrals of monomials of the 
form E"^D{vi)D{v2)---D{vn) and using the relations D{v)\b >= > and < a\E = a < a\. 
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The radius of convergence is also the absolute value of the nearest singularity of / to the origin. We also 
know that all of the coefficients of the Taylor expansion of / in terms of s are positive. This assures that 
the nearest singularity of / lies in fact on positive real half-line. That is, the current at the thermodynamic 
limit is real and positive, as it should be. A similar method works for the average density of particles of each 
species as well. The global density of particles of type i is 

To evaluate this, we use fugacities zq and z{v) to define an operator C[z] as 

C[z] ■.^zoE + J z{v)a{v)D{v)dv. (30) 

Note that we have C[l] = C, where by 1 we mean the fugacities zq = z{v) = 1. It is straightforward to see 
that 

p(.) = fM^ln<a|C-NI^>L^^. 

Once again, the right-hand of this can be expressed is terms of the radius of convergence i?(z) of a formal 
series /(z; s, a, 6) defined as: 

oo 

/(z;s;a,6) ^ < a|C^[z]|6 > . (32) 

N=0 

-1 

Using an equivalent definition for the radius of convergence as R{z) := liniTv >oo ^(< a|C^[z]|6 > ^ , we 

have: 

5 1 

p(w) = z(w)-j---ln— -- . (33) 
OZ(v) n(z) z=l 

So the key step in obtaining the physical quantities is to calculate the functions (|7|) and (|3|) , which we call 
the generating functions for currents and average densities respectively. 
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3.1 Exact calculation of the generating function /(s; a, h) 



First we use ( p3D to obtain a recursion relation for < a|C |& > 



<a|C7^^+^|6> 



^<a\C''\u><-\C\h> 
2mu u 



Au 
27rm 



< a\C'^\u > 



9{b) 



1 



l-h/u 



(34) 



where the function g has been defined in (g4|), i.e. g{h) =< >. Multiplying both sides of ( |34| ) by s and 



summing over N from zero to infinity, we arrive at 



1 f dii 

-[/(s;a,&) -/(0;a,5)] = * —^f{s;a,u) 
s J Imu 



g{b) + - 

u 



l-b/u' 



(35) 



The generating function, which we calculate in this way, will be restricted to the domain \ab\ < 1. After 
calculating it for this region of parameters, we will analytically continue it for other values of parameters 
as well. Since |i| > |u| > |6| the integrand in the right-hand side of ( ^5|) has just two poles inside the 
integration contour; one at tt = and the other at u = b. This is true provided f{s;a,u) itself is analytic 
for u inside the integration contour. However, we know that for small values of its arguments the function 
f{s;a,u) is analytic (see the remark after eq. (p^)). The result of this calculation will be vahd for small 
values of the arguments of the generating function. One can then use analytic continuation to obtain more 
general results. Knowing the non-analytic structure of the integrand, one can use Cauchy's theorem to 
evaluate the right-hand side of (^Sh: 



-[f{s;a,b) - /(0;a,6)] ^ 



f{s;a,b). 



(36) 



Solution of this equation for /(s; a, b) yields 



f{s;a,b) 



sf{s;a,0)-bf{0;a,b) 
b{.s[gib) + l/b] - 1} 



(37) 
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Note that from 

/(0;a,5)=<a|fe>=^^, (38) 

and g{b) = ^^-^-^ is a function which can be determined once the data of the problem (i.e: the distribution 
function P{v) ) are given. Equation ( p7[ ) then suggests that f(s;a,b) is known, provided a two-variable 
restriction of it, namely /(s; a, 0) is known. Eq.(^7|) contains even more information. To see this, notice that 
from (|3^) it seems that there is a pole for s at 

so = Sib) ^ (39) 

3(6) + 1/6 

From the definition of g(b), it is seen that, as b tends to zero, g{b) tends to unity. So, for small values of 6, 
5(6) behaves like 6: 

5(6) - 6 + 0(6^), as 6^0. (40) 

But this means that as 6 tends zero, the radius of convergence for the variable s tends to zero, and this can 
not be the case, due to the remark after (p7|). To avoid this apparent paradox, it must be true that s — 5(6) 
must not really be a pole, at least for small values of 6. This means that the numerator in (|3^) must also 
vanish for s = 5(6). For this to be the case, we should have 

5(6) /[5(6);a,0] = 6/(0; a, 6) = (41) 

This equation allows us to determine the function /(s;a,0) and hence via (|37|), the complete generating 
function. Denoting the inverse function of 5 by -B: 

Sm.s)] = s, (42) 



we find 



f{s;a,0)^- (43) 
si — aB[s) 
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Note, however, that S is not in general one to one and in different domains of its arguments it has different 
inverses. In fact we will show later that S is at most a two-to-onc function with the property S'(O) = 0. By 
the inverse B we mean the one that tends to zero as its argument tends to zero: 



lim B{s) = 0. (44) 

s—>0 



Inserting ( [43| ) in (|37]), and using (41), we find 

B{s) 



l-aB(s) l — ab 
/(s; a, b) = ^ ^ . (45) 



' 1 



_S{b) 

This is the final form of the generating function. For any given probability distribution of hopping rates, 



one can obtain S{b) and hence B{s) from (p4[), and (|39|), which after insertion into ( [45[ ) gives the complete 
generating function. What we will do in the next sections is to carry out an analysis of the singularity 
structure of this function and determine the currents and hence the different phases of our multi-species 
stochastic process. Our results and analysis depend on the general behavior of the functions S{b) and B{s) 
which in turn depend on the distribution of hopping rates. 



4 The total current and the phase structure of the system 
4.1 Properties of the function S 

As it was seen in the previous section, to investigate the properties of the system, one must know the behavior 
of the function S. We have 

^ = i+.W, (46) 

from which we find 
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Thus S~^{b) is a concave function. Combination of this with the fact that this function is positive for 
b e (0,t;i) and tends to infinity as b — > 0, impHes that has at most one minimum in (0,i>i). Hence S 
is a positive function in (0,t;i) and has at most one maximum in this domain. Note also that 5(0) = 0. 
The phase structure depends crucially on whether S attains a local maximum in this domain (i.e.(0, or 
not. This is easily checked from the sign of S'{vi) = (^y-y^yi ^- This quantity is determined only by the 

probability distribution of hopping rates and this is where this function plays its essential role. To emphasize 
the dependence on the distribution function we denote this quantity by l[a] 

lW]--=^-( , ,2 )- (47) 

As we will see if Z[cr] < there are three regions in the phase diagram namely the high-density, the low-density 
and the maximum current phases. On the other hand if l[a\ > 0, the high density phase disappears and only 
the low-density and the maximum current phases remain. We call these two regimes, the 3-phase and 2-phase 
regimes respectively. Qualitatively the transition from the 3-phase to the 2-phase regime is accomplished by 
shifting the distribution function from low speeds to higher speeds. As an example if a{vi) ^ 0( i.e: if there 
is a significant relative probability of injecting slow particles to the system), then it is clear that l[a] = —00, 
which means that we are in the 3-phase regime. The case of discrete values of particle velocities is a special 
case in this category. However if a{v) approaches zero slowly enough as v approaches vi, (i.e: if the chance 
of entrance of slow cars is small), then we will be in the two-phase regime. The exact criterion is given by 
the parameter l[ac\ = 0. 

In the sequel we will need one further property of S. From its definition it is seen that S as a function of the 
complex variable b has a singularity a,tb = vi. If the distribution is discrete, this singularity is a simple pole. 
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Figure 5: The generic form of the function 5(6) that produces the three phase regime 
If the distribution is continuous, S has a branch cut on a segment of the real line beginning from b = Vi to 
+00. To see this we use — tt- = T ii^5{x) to obtain 

4.2 The singularities of the generating function /(s; a, h) 

In order to determine the phases we have to determine the singularities of the generating function. For the 
case l[a\ < where S has a maximum Sm at hm S (0, vi) (Fig. 5), we define two right inverses for S. 

One is the function B defined in the previous section. It is defined in the interval [0,5^], and has the 
following properties. 

' S[B{s)\ =S, 0<S<Sm 

(48) 

^B[S{b)] = b, 0<b<bm 
17 




Figure 6: The functions B{s) and B{s) 

The other function is B, defined in [S{vi), Sm], with the foUowing properties 

S[B{s)] ^ s, S{vi) < S < Sm 

(49) 

B[S{b)] = 6, b,n<b< vi 
The graphs of these functions are shown in figure 6. 

We now consider the generating function as given in (^). The singularities of / as a function of s may 

be of one of following types. The first singularity denoted be Sq. may arise from vanishing of (1 — aB{s)) in 

the numerator. That is 

= - = a, (50) 
a 

or 

Sa = S{a). (51) 

However ( ^0[ ) has a real positive solution for s if and only if a = (1/a) < 6^ (see fig. (6)). Thus is a 
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singularity provided that a < bm- 

The second singularity which we denote hy sp may arise from vanishing of the denominator, i.e: 

Sf3 = S{b) = S{l-P). (52) 

However for b < 6„j, according to ( |4^ ) the numerator vanishes as well, and sp is no longer a singularity, 
which means that sp is a singularity only when b > bm (or /3 < 1 — bm)- 

Finally, B itself becomes singular at s = Sm- Notice that Sm is greater than Sa and sp both, provided these 
latter two exist (see fig. (5)). 

In the next subsection we use these information to determine the phases. 
4.3 The phase structure for / < 

In this case, the generating function has three singular points, namely Sq., S/3 and Sm, each phase (the analytic 
expression of the current) is determined according to which of these singular points are the smallest. 

• The low-density phase (sq < s;5,s„j): 

In this phase which develops when a < bm and S{a) < S{1 — (3), we have: 

J = S'(a), (53) 

• The high-density phase (sp < Sq,s,„): 

In this phase which develops when f3 < 1 — bm and S{1 — (3) > S{a), the total current is given by 

J^S{l-(3), (54) 

and finally 
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The maximum current phase (sm is the only singularity): This phase exists for the rectangular 
area a > bm and /? > 1 — bm- The current is given by 



J 



(55) 



Summarizing we have: 



S{a), 



J={ S{P), 



a < bm, and S{a) < S{1 - (3) 



(3<l-bm, and S{l-P)<S{a) 



(56) 



.Sm = S{bm), a>bm, and P > 1 ~ bm 
The phase diagram is shown in fig. (2). The coexistence curve between the low and high-density phases is 



obtained by the nontrivial solution of the equation S{a) = S{1 — (3). A parametric representation of this 



curve is: 



( a = B{s) 



S{vi) < S < Sn 



(57) 



1/3 = 1- B{s) 

We see that the main features of the 1-ASEP diagram is present. In this regime the multi-species nature of 



the process has only a minor effect. We will see that this exact result is also substantiated by a domain wall 
analysis in accordance with the analysis of B2] . 



4.4 The phase structure for I > 



In this case the generating function has only two singular points, namely Sa and Sm = S{vi). Consequently 
we have only two phases (Fig. (4)). 
• The low density phase in which 



J = S{a), 



(58) 
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and 



• The maximum current phase in which 



J^S{v,). 



(59) 



In summary 



5(a), 



a < vi 



J 



(60) 



a > vi 



The high-density phase has been shrunk and lost, and the injection parameter determines which phase we 
have. This is due to the fact that in this regime the current density diagram is monotonically increasing 
(see section 7) and hence according to the domain wall analysis , only the low density and the maximum 
current phases are expected to exist. Moreover in the maximum current phase everything is controlled by 
the lowest-speed particles. 

The disappearance of the maximum density phase has an interesting implication, for the occurrence of traffic 
jams and its dependence on the bulk parameters beside the boundary ones. From the above analysis one can 
conclude that the maximum density or traffic jam occurs only when there is a critical probability of having 
particles or cars of slow velocities, the exact criteria is given by the parameter I defined above. 

5 Exact calculation of the generating function /(z; s; a, b) and the 
average densities of each species 

Using the same calculation, which led to (^Sf), we obtain 




(61) 
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and from that we arrive at an expression analogous to (p 



/(z;s;a,6) 



zq s /(z; s; a, 0) - b /(z; 0; a, b) 
6{s[5(z;6) + zo/6]-l} 



(62) 



where 



1 - 6/w 



du (t(w) 



1 - 



From (p^), by a reasoning exactly the same as that of section 2, we arrive at 

B{z;s) b 



f{z;s;a,b) = 



1 — a -B(z; s) 1 — a 5 



5(z;fo) 



1 



where 



1 ^0 , ^ 



and i?(z; s) is that right-inverse of S{z; b) which tends to zero as s ^ 0. 



(63) 



(64) 



(65) 



According to (33), or its analogue for the case of continuous distribution, the average density of particles 



of speeds between v and v + dv. 



p{v) — z{v)- — ^— In ■ 



(66) 



5z{v) R{z) z=i 

Knowing the smallest singularity of f{z;s;a,b) is sufficient to obtain the average densities p{v). Once again, 



we can distinguish three phases: the low-density phase, the high-density phase, and the maximum-current 



phase. In the low-density phase i?(z) = S{z,a). Thus p{v) is obtained from (63063) and (B6|) as follows 



p{v) = ziv) 



6z{v) 



-In 



a 



dv cr{v) 



z{v) 



1 — a/v 



a{v)S{a) 
1-^ ■ 



(67) 



The expression for the high-density phase is similar and reads 



p{v) 



<j{v)S{l-f3) 
l-il-P)/v- 



(68) 



The treatment of the maximum current phase however requires more care. In the two-phase regime (i.e: 
I > 0), where the function S does not attain any maximum in [0,111], the maximum current is S'(fi), the 
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situation is the same as above, and we have 



p{v) = 



a{v)S{vi) 
1 — vi/v 



(69) 



In the 3-phase regime (i.e: I < 0) however, the maximum current is S{bm), where bm itself depends on z. 



Therefore, 



+ 



d 1 



b=b„ 



(70) 



dz{v) S[z,bm.{z)] z=i 6z{v) S{z,bm) z=i Sz{v) z.=idbS{b) 
The second term in the right-hand side is, however, zero, since is minimum at bm- So we arrive at the 



expression 



p{v) 



a{v)S{bm) 
1 - bm/v 



(71) 



To summarize, we have 



p{v) 



<7{v)S{x) 

1 — x/v 



(72) 



where 



' a, low-density phase 



1 — /?, high-density phase 



bm, maximum current phase in the three-phase system 



maximum current phase in the two-phase system 



(73) 



The average density of vacant sites po can also be obtained either by using the formula po = zq-^ In 



z=l 



or by using the sum rule 



po+ dv p{v) = 1. 



(74) 



From ([74|), one obtains for each phase determined by the parameter x defined in (|73|) 



Po = 1 - / av- — . = 1 - S{x)g{x) 

\ — x/v 



(75) 



where S{x) = 5(1, x) and g{x) = g{l,x). After using (pq), this gives 



(76) 
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6 Examples 



6.1 The single species ASEP 



In this case a{v) — 5{v — 1). All we need to know to treat this special case is the function S{b) as given in 



eq. 



From (139) we find 



S{b) = 



1 



6(1 - b) 



(77) 



with I ~ Q,b„i = -^j and Sm = j- Thus according to (pfl) we have the following phases in accordance with 



previous results : 



a{l — a) a < f3, and Oi < \ 



a > ^, and P < \ 



a > i and P > h 



(78) 



6.2 Continuous distributions; concrete examples of the disappearance of the 
maximum density phase 

In this section we consider two classes of distribution functions to see concretely the transition between the 
two and three phase regimes. Both of the distributions must be such that they vanish at wi, otherwise as we 
have already remarked I = —oo and we have three phases. . For convenience we also rescale the time so that 
the average velocity is no longer equal to unity. Correspondingly the expression (^^ for I [a] is replaced by 
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where v is the average hopping rate. 

The first distribution that we study has a finite support. 

a{v) = Am{v ~ vi)™-, m > 0, vi<v<V2- 

This kind of distribution has also been considered in a related context by Evans . Note that m need not 
be an integer, and A„ is a normalization constant to be determined shortly. It is convenient to evalute the 
following integrals: 

oo k < —1 

4:= / (t. -t.,)*<ii.= { . (80) 



Obviously A,n — j—- Simple calculations also yeild: 



w = f 1 H (81) 

m + 2 



and 



= (-)' - -^{Im-l + ViIm-2), (82) 



from which we find 



l[a] = { ^ ^ . (83) 

1 < m 



-oo < m < 1 

l(2ii,m,) = (^)2 - (to + l)v( 1 . + . ^^fi ^) 

The sign of the parameter l{vi, V2, m) determines if the maximum density phase exists or not. In an extreme 
case the analysis of this quantity is quite simple. For very large to, we find after inserting (Is^) in (^) 



l{x := — , to) = x^{1~ ^j^) - —x{x - 1) + 0{m-^). (84) 
vi X ^ 1 m 

which implies that to zeroth order, if a; > 2 (i.e. V2 > 2wi), there is no maximum-density-phase. To first 

order the above condition is modified to V2 > '2vi{l + ^). 

We now consider another distribution with infinite support. 
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cr(-y) = - ui)™e to > 0, Vi < v 



The analysis is similar to the previous case. We have 



A: < -1 



from which wc find 



or 



Jk / (v^vife —dv^ { (85) 

Jfc A'=+ir(fc + 1) -Kk 



V — vi + toA (86) 



= ( — )' - -J-{Jm-l + ViJm~2). (87) 



— oo < m < 1 

l(x:=^,m) Km 



where 



l(., TO) = x(2 + x) - 4 ( "^+(^"\-^)" ) . (89) 

\ TO — 1 



. Again to zeroth order of — , we find 

o rn, ' 



a; a; 



/= (.T--)(a; + - + 2) (90) 



which implies that when ^ > 1, the high-density-phase disappears. 



6.3 The p-species ASEP, with a hopping rate much lower than the others: 

The case of a fixed or moving impurity has been studied in many previous works as for example in [ p9| pO| , |3l[ . 
In the present framework we can consider a new case where the number of impurities is not one or even 
fixed. That is we allow very slow particles to have a chance of entering into and leaving the system. That 
is we take 

1 ^ 
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High density phase 



1-v, 



a 



b 



m 



Figure 7: Phase diagram of the p-species ASEP when one of the hopping rates is much smaller than the 
others. 

and let one of the particles has a speed much lower than the rest: that is : b < vi «< V2 < V3 < ■ ■ ■ < Vp. 
Then with the approximation 1 — ^ ~ 1, for i = 2,3,- ■ -p, one can write: 



i=l Vi ' VI 



(91) 



from which one obtains: 



S-\b)^-+g{b) = 



6 + 1 



b p{vi — 6) 



(92) 



For this function we have: 



bm = 



Vl 



and s„ 



Vl 



1 + 



l + vi+2jf 



(93) 



Besides vi , a and /?, only the number of species plays a role here. 
The phase diagram is shown in Fig. (7). 

The following features are readily observed. Compared to the other two phases, the size of the low density 
region is very small, as we expect on physical ground. Only for very small injection rate a and for very large 
extraction rate f3 can this phase exist in the system. Even in the maximal current phase the current which is 
given by Sm is seen to be small and limited by the speed of the lowest particles. For fixed a, f3, and vi, 
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as the number of species p increases, the value of bm approaches vi and hence the high density phase begins 
to shrink, leaving only two phases in the system. 

It is instructive to calculate the relative numbers of particles of different types including holes in each of the 
above phases. From (^2|) we find 

pM. - ^"^ (94) 

where in each phase x is given as in ([73|). Inserting the relevant values of x and in the approximation 
-^^ « 1 Vi > 2, we find: 



1_ a. ■ 



Low-density phase 



1 - -f , High-density phase (95) 



1 — maximal current phase 
In the low density region where a takes values from to this ratio takes values from 1 to 1 



vi 



1 Thus in this region almost all types of particles are present in the system. In the high density 



region this ratio is at most equal to 1 which is determined by the interplay of the lowest speed 

and the number of species. For large p, it is indeed a very small value , indicating that the system is almost 
filled by the lowest particles 

This ratio takes its maximum value in the entire maximal current phase. One can also obtain the ratio 
of the density of lowest speed particles to that of the holes. From (^2|) and ( |76| ) one obtains 

1 a 

1 b , . 

pT^ (96) 



pM 

Po 



1 b„ 



p i_ 



where from top to bottom we have listed the low density, the high density and the maximal current phases. 
Again we can find the limiting values of this ratio, to see how crowded the system is, in each phase. It is 
simple to see that the ratio ranges from to in the low density, and from to infinity, in the high 

28 



density and is fixed at y in the maximal current phase. 

With a httle more effort, starting from (|9^), the coexistence line between the low-density and the high-density 
phases is found to be given by: 

P = a + l- Pf+''y'-P^ (97) 
[a + pjvi — ap 

7 Discussion 

The results that we have obtained on the phases and currents are exact. We can get a feeling for these 



results, based on the intuitive arguments of domain wall dynamics [^, 33, g9|, 16, What we will do in 
this section is to formally adopt the analysis of and redrive our exact results. The essential result of 
is that for all single species processes which have single peak current density relation, the phase diagram 
of the ASEP is generic, that is the possible phases are the low density, the high density and the maximal 
current phases. Roughly speaking one expects that for a small and (3 large, the low density phase denoted 
schematically by (000000000), prevails in the system, and for a large and /3 small, the high density phase 
denoted by (11111111) prevails. However when there is no restriction on the injection and extraction rates 
of the particles, that is for a and (3 large, the current reaches its maximum value allowed in the current 
density diagram, this new phase being called the maximal current phases and denoted by (mmmmmmmm). 
The exact shape of the phase diagram and the coexistence lines are obtained by studying the dynamics of a 
supposedly formed domain wall at sufficiently late times between any pair of these phases under appropriate 
conditions. For example when a is small and /3 is large, the late time configuration is supposed to be 
(000000111111). It is also assumed that deep into each of the two segments we have a product measure 
with constant density. This assumption is well-founded M, RSl E9, by numerical, mean field and exact 
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solutions . The velocity of such a domain wall is then given by the formula 

y^Jo_^^ (98) 

Po - Pi 

where Jq and po ( resp. Ji and pi ) are the current and density to the far left ( resp. right ) of the 
domain wall. The sign of this velocity determines the prevailing phase and setting this velocity equal to zero 
determines the coexistence line. In the latter case the two phases coexist due to dominance of fluctuations 
in the rms position of the domain wall. For the currents and the densities in (|9^ ) one uses the mean field 
values. For the maximal phase, one uses the density which maximizes the current in the current density 
diagram. 

The above analysis can be readily applied to the multi-species case. On the assumption that the coarse 
grained bulk current is given by the uncorrelated Bernoulli measure [^^, we can use the one-parameter 
family of one dimensional representations for the bulk relations in (|[^) to obtain |Q 

D{v) = ^ (99) 
V — b 

and consequently the following forms for current and total density: 



Note that the right hand side of ( |100| ) is exactly the function 5(6) defined in equation (|39|). However, before 
using equation (p8|), we need to determine the free parameter b and its range, in the Bernoulli measure for 
each phase. For the maximum current phase the parameter is obviously bm which maximizes J{b). This 
is exactly the parameter, which has been defined in section (4.1). For the other two boundary-controlled 
phases, the parameter b should be fixed by compatibility with the boundary conditions of (|§-^, according 
to which a product measure coupled to a left reservoir injecting particles at rate a should have b — a and a 
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product measure coupled to a right reservoir extracting particles at rate (3 should have 6=1-/3. Instead 
of using this type of argument which is based on MPA relations one can follow the more general argument 
suggested in ||l^, to match the boundary rates with the bulk densities. 

Note also that for the low density phase b < bm and for the high density phase b > Thus we have 
Jo = Jib — ct), Ji = J{b =1-/3) and Jm = J{bm)- Equating the currents we obtain exactly the phase 
structure previously obtained by exact solution. Moreover the size of the maximum density region in the 
phase diagram depends on the value vi — b^ (see Fig. (5)). When vi — b,n approaches zero, the size of this 
region shrinks and we remain only with two phases. This is again in accord with our exact solution. 
To conform completely to the picture advocated in |^ we should have discussed various phases according 
to the behavior of the function J{p) and not the function J(b). However in our case the qualitative behavior 



of these two functions resembles each other. In fact it is seen from (IOC -101) that p is a monotonically 
increasing function of 6, which attains its maximum pi at b — vi. Moreover J{p = 0) = 0, J'{p = 0) = 1, 
and finally | |23[ | J is a convex function of p. To see if J(p) attains a local maximum in its domain of definition 
[0, pi] or not we evaluate ^ at p — pi and find 

1 



= ¥^ - ^ 7, (102) 



1 



\ v — v-i 

Thus here also the value of the parameter / determines the answer to the above question. 
We should stress that the above arguments due to their qualitative nature Is not by no means a substitute 
for exact solutions. However it is remarkable that in view of the crude approximations involved, they can 
predict exact results. 
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9 Appendix 

For our presentation to be consistent, we have to show how the infinite dimensional algebra is derived 
in the MPA formalism. This can be simply done by a slight modification of the relations in j2^. In the 
general case the Hilbert space of each site of the lattice which we denote by h is generated by a discrete 
state |0) (when the site is empty) and a continuous set of states |u), u £ (0, oo) (when the site is occupied by 
a particle of intrinsic velocity v). We denote these states by different symbols to avoid confusion with the 
states of the representations of the algebra. The states are normalized as : 



(0|0) = 1 



(0|u) = {v\0) = 



{v\v') ^ S{v -v') 



(103) 



The Hamiltonian is 



k=N-l 




h 



B 

'k,k+l 



(104) 



where h 



is given by 



h 



B 




(105) 



The boundary Hamiltonians hi and h]\f are: 




(106) 




(107) 
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Note that the distribution a{v) only enters hi, which points to the fact that the distribution a{v) refers to 
the particles injected to the system. 

Inserting these Hamiltonians in the standard formulas of the MPA, i.e: 

h^A(^A = X (giA-A(g>X <W\hiA + X = HnA- X\V >=0 (108) 

with the following form of the auxiliary vectors A and X 

A = E\0) + J cr{v)D{v)\v)dv 

X = -\0) + J <7{v)v\v)dv (109) 

leads to the algebraic relations (§-||). Note that A and X are operator valued vectors in the Hilbert space 
of one site of the lattice, as they should be in the MPA formalism. 
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Figure Captions 

Fig. 1: Phase diagram for the single species ASEP. 

Fig. 2: Phase diagram for the multi-species ASEP, when l[a] < 0. 

Fig. 3: Two phase diagrams for the multi-species ASEP for different distributions of hopping rates. 

In both cases l[a] < 0. 

Fig. 4: Phase diagram for the multi-species ASEP when l[a] > 0. 

Fig. 5: The generic form of the function S(b) that produces the three phase regime. 

Fig. 6: The functions B{s) and B{s) . 

Fig. 7: Phase diagram for the p-species ASEP, when one of the hopping rates is much smaller than 
the others. 



37 



